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Abstract. The aim of this paper is to give a higher dimensional equiva- 
lent of the classical modular polynomials {X, Y) . If j is the j-invariant 
associated to an elliptic curve Ek over a field k then the roots of $e{j-, X) 
correspond to the j-invariants of the curves which are ^-isogeneous to 
Eh. Denote by Xo{N) the modular curve which parametrizes the set of 
elliptic curves together with a A'^-torsion subgroup. It is possible to inter- 
pret ^t(X, Y) as an equation cutting out the image of a certain modular 
correspondence Xo{l) — > A'o(l) x A"o(l) in the product A'o(l) x A'o(l). 
Let 3 be a positive integer and n £W . We are interested in the moduli 
space that we denote by Mn of abelian varieties of dimension g over a 
field k together with an ample symmetric line bundle ^ and a symmetric 
theta structure of type n. If £ is a prime and let Z = [t, . . . ,tj, there 
exists a modular correspondence Mj^ — > M— x M—- We give a system 
of algebraic equations defining the image of this modular correspondence. 
We describe an algorithm to solve this system of algebraic equations 
which is much more efficient than a general purpose Grobner basis algo- 
rithm. As an application, we explain how this algorithm can be used to 
speed up the initialisation phase of a point counting algorithm. 
Keywords: Abelian varieties, Theta functions, Isogenies, Modular cor- 
respondences. 



1 Introduction 

The aini of this paper is to give a higher dimensional equivalent of the classical 
modular polynomials (l>i{X,Y). We recall that 'I>i{X,Y) is a polynomial with 
integer coefficients and that if j is the j-invariant associated to an elliptic curve 
Efe over a field k then the roots of <l>£{j,X) correspond to the j-invariants of 



elliptic curves which are ^-isogeneous to Ek- These modular polynomials have 
important algorithmic applications. For instance, Atkin and Elkies (see |Elk98| ) 
take advantage of the modular parametrisation of ^-torsion subgroups of an 
elliptic curve to improve the original point counting algorithm of Schoof [Sch95| . 

In [SatOOj , Satoh has introduced an algorithm to count the number of rational 
points of an elliptic curve Ek defined over a finite field k of small characteris- 
tic p which rely on the computation of the canonical lift of the j-invariant of 
Ek- Here again it is possible to improve the original lifting algorithm of Satoh 
|VPV01|LL06| by solving over the p-adics an equations given by the modular 
polynomial $p{X,Y). 

This last algorithm has been improved by Kohel in [Koh03| using the notion 
of oriented modular correspondence. For N E N* , the modular curve Xo{N) 
parametrizes the set of isomorphism classes of elliptic curves together with a 
TV-torsion subgroup. For instance, the curve Xq{1) is just the line of j-invariants. 
Let p be prime to N. A rational map of curves Xd{pN) Xo{N) x Xq{N) is an 
oriented modular correspondence if the image of each point represented by a pair 
{E, G) where G is a subgroup of order pN of i? is a couple {{Ei, Gi), {E2, G2)) 
with El — E and Gi is the unique subgroup of index p of G, and E2 = E/H 
where H is the unique subgroup of order p of G. In the case that the curve, 
Xq{N) has genus zero, the correspondence can be expressed as a binary equation 
y) = in Xo{N) X Xii{N) cutting out a curve isomorphic to Xo[pN) 
inside the product. For instance, if one consider the oriented correspondence 
-'^o(^) -'^o(l) X -''^0(1) for I a prime number then the polynomial defining its 
image in the product is the modular polynomial '^eXX, Y). 

In this paper, we are interested in the computation of an analog of oriented 
modular correspondences for higher dimensional abelian varieties over a field k. 
We use a moduli space which is different from the one of [Koh03| . We fix an 
integer g > for the rest of the paper. In the following if n is an integer, n 
denotes the element (n, . . . ,rt) S . We consider the set of triples of the form 
{AkT^, On) where Ak is a g dimensional abelian variety equipped with a sym- 
metric ample line bundle .jSf and a symmetric theta structure On of type n. Such 
a triple is called an abelian variety with a ri-marking. To a triple [Ak,■^,On)^ 
one can associate following |Mum66| its theta null point. The locus of theta null 
points corresponding to the set of abelian varieties with a n-marking is a quasi- 
projective variety Mn- Moreover, it is proved in ^Mum67j that if 8|n then Mn 
is a classifying space for abelian varieties with a n-marking. We would like to 
compute oriented modular correspondences in AAn- 

For this, let {Ak,^ ,0j:^) be an abelian variety with a (£n)-marking. We 
suppose that ^ and n are relatively prime. From the theta structure Oj^, we 
deduce a decomposition of the kernel of the polarization K{^) = Ki{^) x 
K2{^) into maximal isotropic subspaces for the commutator pairing associated 
to ^. Let K{^)[e] = Ki{^)[£] x K2i^)[£] be the induced decomposition of the 
^-torsion part of K{Ji'). Let Bk be the quotient of Ak by K2{^)[(\ and Ck be the 
quotient of Ak by Ki{^)[£]. In this paper, we show that the theta structure of 
type in of Ak induces in a natural manner theta structures of type n on Bk and 



Cfc. As a consequence, we obtain a modular correspondence, Mj^ — > Mn x Mn- 
In the projective coordinate system provided by theta constants, we give a system 
of equations for the image of TM-^ in the product Ain x as well as an efficient 
algorithm to solve this system. 

This paper is organized as follows. In Section[2]we recall some basic definitions 
and properties about algebraic theta functions. In Section |31 we define formally 
the modular correspondence, and then in Section|3]we give explicit equations for 
the computation of this correspondence. In particular, we define a polynomial 
system (the equations of the image of A^t;^), which solutions give theta null 
points of isogeneous varieties. In Section [5l we describe the geometry of these 
solutions. The last Section is devoted to the description of a fast algorithm 
compute the solutions. 



2 Some notations and basic facts 



In this section, we fix some notations for the rest of the paper and recall well 
known results on abelian varieties and theta structures. 

Let Ak be a (7 dimensional abelian variety over a field k. Let ^ be a degree 
d ample symmetric line bundle on Ak- From here, we suppose that d is prime 
to the characteristic of k or that Ak is ordinary. Denote by K{^) the kernel 
of the polarization and by G{^) the theta group (see [Mum66| ) associated 
to ^ . The theta group G{J^) is by definition the set of pairs (a;, ip) where x 
is a geometric point of K{^) and if) is an isomorphism of line bundles ip : 
^ T*I£ together with the composition law (a;,-)/') o (y, <p) = (x + y,T*^p o 
(f). Let S = (di, . . . ,dg) be a finite sequence of integers such that di\di^i, we 
consider the finite group scheme Z{6) — {'Z,/diZ)k Xfc ... x^ {'Z/dgZ)k with 
elementary divisors given by S. For a well chosen unique S, the finite group 
scheme K{6) = Z{6) x Z{6) (where Z{d) is the Cartier dual of Z{S)) is isomorphic 
to i^(-Sf) (see IMumTOj '). The Heisenberg group of type 6 is the scheme H^S) = 
Gm,k X Z{S) X Z{S) together with the group law defined on geometric points by 
{a,xi,X2)-iP,yi,y2) = {a.(3.y2ixi):Xi +yi:X2 +2/2)- We recafi fMum66" that a 
theta structure 0s of type 5 is an isomorphism of central extension from Ti.{d) 
to G{^) fitting in the following diagram: 



0- 



m.k ■ 



Kid) 







(1) 



06 



m.k ' 



G(^) 



if(if) 



We note that 0s induces an isomorphism, denoted 0s in the preceding di- 
agram, from K{6) into K{^) and as a consequence a decomposition K{^) = 
Ki{^) X K2{^) where K2{^) is the Cartier dual of Ki{^). The data of a triple 
(^fe, 0s) defines a basis of global sections of .Sf that we denote (i?j)jez(5) ^^id 
as a consequence an morphism of Ak into P^T^ where d = 0^=1 ^-i is the degree 
of We briefiy recall the construction of this basis. We recall |Mum66t pp. 



291] that a level subgroup K of G(^) is a subgroup such that K is isomorphic 
to its image by k, in K{^) where k is defined in ([1]). We define the maximal 
level subgroups Ki over Ki {^) and K2 over K2 {^) as the image by 0s of the 
subgroups [ItX,Q)x^z{5) and (1,0,2/)^^^^^^) of T-L[5). Let be the quotient of 
Ak by 7^2 (-S^') and n : Ak ~* A^ be the natural projection. By the descent the- 
ory of Grothendieck, the data of K2 is equivalent to the data of a couple (^0, ^) 
where is a degree one ample line bundle on A^, and A is an isomorphism 
A : 7r*(.ifo) I£ . Let sq be the unique global section of S£q up to a constant 
factor and let s = A(7r*(so)). We have the following proposition (see [Mum66| l 

Proposition 1. For all i e Z{S), let {x^,ip^) ^ ©^((l, z, 0)). We set df' = 

(Tl^.tpiis)). The elements {'&f^)i(zz{S) form a basis of the global sections of 
which is uniquely determined up to a multiplication by a factor independent of i 
by the data of 0s . 

If no ambiguity is possible, we let i?f = for i G Z{5). 

The image of the zero point of A^ by the morphism provided by 0s , which 
has homogeneous coordinates {■di{0))i^z{5)^ is by definition the theta null point 
associated to {Ak,^,0s)- If 0s is symmetric |Mum661 pp. 317], we say that 
{Ak,.Si',0s) is an abelian variety with a 5-marking. The locus of the theta null 
points associated to abelian varieties with a (5-marking is a quasi-projective va- 
riety denoted Ms- 

Let {Ak,^,0s) be an abelian variety with a J-marking. We recall that the 
natural action of G{^) on the global sections of ^ is given by {x, ■;/')■/ = ^Ixi^if) 
for / £ r{J^) and (a;, %p) e G{^). There is an action of 1-L{5) on {'&i)iez(S) given 
by: 

{a,i,j).{>k = aes{k + i,-j)-dk+i, (2) 

for {a,i,j) S Ti-^S) and es the commutator pairing on K(S), which is compatible 
via 0s with the natural action of G(^) on {'&i)i(rz{S) ■ Using one can compute 
the coordinates in the projective system given by the {Oi)i^z(S) of any point of 
K{^) from the data of the theta null point associated to (A^, 0s)- 

Let 5 = {61, ...,5g) e W and 5' = {5[, ...,5'g) € W, S\d' means that for 
i = 1,. . . ,g, Si\S-. If n e N, n\6 means that (n, . . . , n) e W\S. If 6\S' we have the 
usual embedding 

i : Z{5) Z{5'), (a;i)i£{i,...,g} 1-^ {S'JSi.Xi) (3) 

A basic ingredient of our algorithm is given by the Riemann relations which 
are algebraic relations satisfied by the theta null values if A\5. 

Theorem 1. Denote by Z(2) the dual group of Z{2). Let {ai)i,^z(S) ^6 theta 
null points associated to an abelian variety with a 5-marking [A^,^ ,0s) where 
2\5. For all x, y,u,v £ Z{26) which are congruent modulo Z(2), and all x G Z{2), 



we have 




Here we embed Z{2) into Z{6) and Z{S) into Z{26) using ([3]). 

It is moreover proved in |Mum66| that if A\S the image of by the projec- 
tive morphism defined by Os is the closed subvariety of V^"^ defined by the 
homogeneous ideal generated by the relations of Theorem [TJ 

A consequence of Theorem [T] is the fact that if 4|5, from the knowledge of a 
valid theta null point {ai)i(zz{5), one can recover a couple (A^,^) from which 
it comes from. In fact, the abelian variety Ak is defined by the homogeneous 
equations of Theorem[TJ Moreover, from the knowledge of the projective embed- 
ding of Ak, one recover immediately ^ by pulling back the sheaf 0(1) of the 
projective space. 

An immediate consequence of the preceding theorem is the 

Theorem 2. Let {ai)i^z{S} be the theta null point associated to an abelian vari- 
ety with a 6-marking [Ak, ^ , Os) where 2\S. For all x,y,u,v £ Z(2S) which are 
congruent modulo Z{S), and all x G Z{2), we have 



As Os is symmetric, the theta constants also satisfy the additional symmetry 
relations ai ~ a^i, i G Z{S). 

The Theorem [5] gives equations satisfied by the theta null points of abelian 
varieties together with a (5-marking. Let Ais be the projective variety over k 
defined by the symmetry relations together with the relations from theorem [21 
Mumford proved in |Mum67] the following 

Theorem 3. Suppose that 8\S. Then 

1. M.S is a classifying space for abelian varieties with a 5-marking: to a theta 
null point corresponds a unique triple {Ak,^,Os)- 

2. Ms is o,n open subset of Ms- 

A geometric point P of Ms is called a theta constant. If a theta constant P 
is in 5 we say that P is a valid theta null point, otherwise we say that P is a 
degenerate theta null point. 

Remark 1. As the results of Section [5] show, Ms may not be a projective closure 
of Ms- Nonetheless, every degenerate theta null point can be obtained from a 
valid theta null point by a "degenerate" group action (see the discussion after 
Proposition [7]), hence the notation. 




te2(2) 




tez{2) 



tez(2) 



3 Theta null points and isogenics 



In this section, we are interested in the following situation. Let £ and n be 
relatively prime integers and suppose that n is divisible by 2. Let (A^, Jf, 0^;^) 
be a f;-dimensional abelian variety together with a (£n)-niarking. We recall that 
the theta structure induces a decomposition of the kernel of the polarization 

K{^) = i^i(^) X K2{^) (4) 

into maximal isotropic subgroups for the commutator pairing associated to ^ . 
Let K he a. maximal isotropic ^-torsion subgroup of K{^) compatible with 
the decomposition (j4|). There are two possible choices for one contained in 
Ki{^), the other one in K2{^)- In the next paragraph, we explain that a choice 
of K determines a certain abelian variety together with a n-marking. The main 
results of this Section are Corollary [1] and Proposition [3] which explain how to 
compute the theta null points associated to the abelian variety together with a 
n-marking defined by a choice of K. 

Let Xk be the quotient of by K and let tt : Afc — > Xk be the natural 
projection. Let k : G{^) K{^) be the natural projection deduced from the 
diagram ([1]). As ivT is a subgroup of K{^), we can consider the subgroup G of 
G{^) defined as G = k^^{K). Let K be the level subgroup of G(^) defined as 
the intersection of G with the image of (1, x, y)^x,y)ez(hi)xz(i^) *^ 'Hi^n) by 0-^. 
By the descent theory of Grothendieck, we know that the data of K is equivalent 
to the data of a line bundle ^ on Xk and an isomorphism A : n*{^) J^f. 

Now, we explain that the (£n)-marking on induces a n-marking on Xk- 
Let G*(^) be the centrahzcr of K in G(if). Applying |Mum66[ Proposition 2 
pp. 291] , we obtain an isomorphism 

G*i^)/K ~G{.^) (5) 

and as a consequence a natural projection q : G*{^) G{^'). 

As Tiiji) is generated by the subgroups 1g„ x Z{n) x Qz(n) ^^"^ x ^z{n) x 
Z(n), in order to define a theta structure 0^ : Tiin) G(J?^), it is enough to 
give morphisms 1g,„ x Z{n) x O^^--^ G( JT) and 1g„ x Oz{n) x Z{n) G( JT). 
Let Z*(in) be such that 1g„ x Z*(hi) x O^^j;^) = 0:^{G*{^)) n Z(hi) and let 

Z*(e^) be such that 1g,„ x O^ihl) x z*(e^) = e^{G*i^)) n z(e^). 

As Z*{£n) is in the orthogonal of 0j^{K) for the commutator pairing, we 
have Z*{(.n) = Z{ln) or Z*{ln) = Z{n) depending on the choice of K. In any 
case, there exists a natural projection p : Z*{in) — > Z(n). In the same way, 
Z*{£n) = Z{£n) or Z*(in) — Z(n) and there is a natural injection i : Z{n) 
Z*(E). 



We can define On as the unique theta structure for ^ such that the foUowing 
diagrams are commutative 



(6) 



(1,^,0)^^^.(5^) G*(^) , (7) 



where i is deduced from i and p is deduced from p. Using the fact that is 
symmetric, it is easy to see that On is also symmetric. 

We say that the theta structures and 0n are 7r-compatible (or compati- 
ble) if the diagrams ^ and ([7]) commute. 

Let Ki and K2 be the maximal ^-torsion subgroups of respectively Ki{^) 
and K2{^)- By taking K = K2 and K = Ki in the preceding construction, we 
obtain respectively {Bk,^f),&n) and {Ck,^,0^) two abehan varieties with a 
n-marking. As a consequence, we have a well defined modular correspondence 

■■ -^MnX Mn- (8) 

Let TT : Ak ^ Bk and tt' : Ak ^ Ck be the isogenics deduced from the 
construction. Let [£] be the isogeny of multiplication by £ on Bk and let tt : 
Bk Ak be the isogeny such that [£] = tt o tt. From the symmetry of .if 
we deduce that is symmetric and by applying the formula of [Mum66[ pp. 
289], we have [£]*^q = . The following dia gram shows that Ck is obtain by 
quotienting Bk by a maximal isotropic subgroup of {Bk,^^ ) of order P^. 



Bk . (9) 




Bk Ck 

The following two propositions explain the relation between the theta null 
point of {AkjJif, Oj^) and the theta null points of .ifo. On) and [Ck,^ , O^). 
Keeping the notations of the previous paragraph, we have 

Proposition 2. Let [Ak^^ ,0-^), {Bk,^o,On) and tt : Ak ^ Bk be defined as 
above. There exists a constant factor to G k such that for all i G Z(Ti), we have 



Tr*{'df") — Lwdf In this last relation, Z{n) is identified as a subgroup of Z{£n) 
via the map x ^-^ Hx. 

Proof. This proposition is a particular case of the isogeny theorem [Muni66 . Th. 
4] but we give here a direct proof. 

Let Ka be the level subgroup of G(.2') defined by the image of (1, 0, y)y^z(en) 
by Oj:^ and let Ka be the subgroup of Ak which is the image of iO,y)y^z(in) 
by Oj^. Let Dk be the quotient of Ak by Ka and tta ■ Ak Dk the natural 
projection. The data of Ka gives a couple (.5fA,Ayi) where ^a is a degree one 
line bundle on Dk and \a is an isomorphism Xa ■ t^*a{S£a) We recall that 

K be the level subgroup of G'(.Sf) defined as the intersection of G — k^^(A') 
with the image of {^■,x,y)(^^^y^^z^.^^z(hi) C W(^) by Oj;^. 

In the same manner, we can consider Kb the level subgroup of G(^o) defined 
by the image of (Ij 0, y)^^^^^-— -j by and Kb the subgroup of Bk which is the 

image of ^,v)y(zz{n\ by ©n- By © Kb = t^{Ka) and by definition of tt its 
kernel K is contained in Ka- We deduce that Dk is the quotient of Bk by Kb and 
•ka = ttb OTT where tt^ is the natural projection Bk Dk- Because of (|6]) and the 
fact that Z*{£n) — Z{£n), we have an isomorphism Kb — Kb/K and the data 
of i^s gives a couple {J^b^Xb) where As is an isomorphism As : ■k*b{J^b) — * 
and we have J^b = and \a°t^*a = X o -k* o \b o -k*^. 

If So is the unique global section of I£a up to multiplication by a constant fac- 
tor, we have \a{i^*a{sq)) = A(7r*(Ai3(7r^(so)))). By definition, 'd^ = Xb^i^b^^q)) 
and — Xa(j^*a{,sq))- As a consequence, there exists uj G k such that we have 
that = ujd^~. 

Let s = and s' = t?^'^. We set for all i e Z(£n), {x^, = 6)j;;-((l, i, 0)) 
and for all i € Z(n), (x^, V-) = e^((l,i,0)). Then ^*(i9f'^) = 7r*(7/;^T:^,(s')) = 
^iTl^.TT* {s') by the commutativity of ([7]). But we already know that 7r*(s') = 
and ipiTl^. (ujs) = LO'df'^" . This concludes the proof. 

As an immediate consequence of the preceding proposition, we have 

Corollary 1. Let (A^,^, Oj^) and {Bk,.S^o, On) be defined as above- Let {o^u)u^z(J7i) 
and (bu)uez(-n) be theta null points respectively associated to {Ak, ,Oj^) and 
(i?A;,^o, G'n)- Considering Z{n) as a subgroup of Z{£n) via the map x i— > Ix, 
there exists a constant factor uj €z k such that for all u G Zin), fe„ = ^Ou- 

Proposition 3. Let {Ak,.^,Oj:^) and (G^, .5fo, G^n) be defined as above- Let 

{^u)u^z(In) ^"■^ (c«)tiez(ri) be the theta null points respectively associated to 
{Ak,^, Oj^) and {Ck,^o, On)- We have for all u e Z{n), 

c« = ^ Ou+t, (10) 
tez(e} 

where Z(n) and Z{£) are considered as subgroups of Z(£n) via the maps j i— > £j 
and j nj . 



Proof. The theta structure O-pp^ (resp. 0!^) induces a decomposition of the kernel 
of the polarization X(^) = Ki{^) x K2{^) (resp. K{^) = Ki{^) x K2{^)). 
Denote by K' the kernel of tt'. We have that K' is a subvariety of _ft'i(^) and 
we have an isomorphism: 

a : Ki{.^)/K' ^ Ki{^). 

The hypothesis of [Mum66i Th. 4] are then verified and Equation pO|) is an 
immediate application of this theorem. 



4 The image of the modular correspondence 

In this section, we use the results of the previous section in order to give equations 
for the image of the modular correspondence . 

We let {Bk,.SfQ,0n) be an abelian variety together with a n-marking and 
denote by {bu)uez(n) its associated theta null point. Let be the 2-adic valua- 
tion of n. Unless specified, we shall assume that ^ 3. Let be the reduced 
subvariety of Mn x Mn which is the image of <l>i (Mj;^) in Mn x Mn given on 
geometric points by tt : {au)^f,z{W^) ^ ((a«)uez(n), (Etez(Z) "'u+t)uez{n))- 

Denote by pi (resp. P2) the restriction to 'rf of the first (resp. second) projec- 
tion from Ain X Mn into A4ni and let tti = pio tt, it2 = P2° tt. We would like to 
compute the algebraic set t^2{'^i'^ {{bu)u<az{n)))- We remark that this question is 
the analog in our situation to the computation of the solutions of the equation 
X) defined from the modular polynomial and j G fc a certain j-invariant. 

Let p^(^") = Vvo]{k[xu\u e Z{in)]) be the ambient projective space of M-^^ 
and let / be the homogeneous ideal defining J^j^ which is spanned by the 
relations of Theorem [21 together with the symmetry relations. Let J be the 
image of / under the specialization map 



6„ , if u G Z{n) 



fc[xtj|u G Z{£n)] — !■ k[xu\u G Z{hi), nu 7^ 0], Xu <- 

and let Vj be the afRne variety defined by J. 

Let n° : Pf^"") ^ Pf^"^ and : Pf^""^ ^ Pf be the morphisms 
of the ambient projective spaces respectively defined on geometric points by 
Muez(i^) ^ (a«)«G2(«) and {au)u£Z{I^) ^ (Etez(I) a«+t)«GZ(«)- Clearly, tti 
and 1T2 are the restrictions of 7r° and to Mj;^. The morphism ttj' (resp 772) 
restricts to a morphism tti : Mj;^ Mn (resp tt2 : M-^ Mn) . By definition 
of J, we have Vj = Tri^{bu)uGZ{n)- 

Let S = k[yu, Xv\u G Z{n), v G Z{£n)], we can consider J as a subset of S via 
the natural inclusion of fc[xti|u G Z{£n)] into S. Let C be the ideal of 5 generated 
by J together with the elements j/„ — J2tez(e) ^u+t and C = L' C\ k[yu\u G Z{n)]. 
Let 1/c be the subvariety of A'^'"' defined by the ideal £. By the definition of £, 

is the image by tt2 of the fiber Vj, so that Vc = 7T'2(7rr^(^u)ue2(n))- 



Proposition 4. Keeping the notations from above, let {bu)u£Z{n} be the geo- 
metric point of A4n corresponding to {Bk,^o, On) ■ The algebraic variety — 
Tr2i'^i^ {bu)uez{n)) has dimension and is isomorphic to a subvariety ofVc- 

Proof. From the preceding discussion the only thing left to prove is that has 
dimension 0. But this follows from the fact that the algebraic variety Vj has 
dimension |CL08| . 

From an algorithmic point of view, the hard part of this modular correspon- 
dence is the computation of Vj ~ T^i^{{bu)uez{n)), the set of points in Vj that 
are valid theta null points. We proceed in two steps. First we compute the solu- 
tions in Vj using a specialized Grobner basis algorithm (Section 16. 3p and then 
we detect the valid theta null points using the results of the next section (see 
Theorem |4]). But at first we recall the geometric nature of Vj given by Section [3j 

Proposition 5. Vj is the locus of theta null points {o-u)uez(In) smc/i 
that if {Ak,^,Oj^) is the corresponding variety with a (in) -marking then Oj^ 
is compatible with the theta structure On of B^. 

Proof. Let {o-u)^^z(In) ^ geometric point of Vj. Let (A^,^, O-p^) be a corre- 
sponding variety with (i'n)-marking. If we apply the construction of Section [3l 
we get an abelian variety {B'^,J^q,0^) with a n-marking and an isogeny tt : 
Ak B'f. such that Oj^ is compatible with O^. By definition of J, Corollary [1] 
shows that the theta null point of B' is {bu)u£Z{n)- As > 2, by Proposition 
[T]this means that B' ^ B. Since z/ ^ 3, we even know by Theorem [3] that the 
triples (B^j^o'i ^tt) ^^'^ iBk,^Q, On) are isomorphic, so that Oj;^ is compatible 
with On. 

We say that the isogeny from Section ^Ak-^Bk is the ?-isogeny associated 
to the (^n)-marking of A^. 

5 The solutions of the system 

This section is devoted to the study of the geometric points of Vj. Our aim is 
twofolds. First we need a way to identify degenerate theta null points in Vj, 
and then we would like to know when two geometric points in Vj correspond to 
isomorphic varieties. 

If {ciu)y^^z(hi) ^ valid theta null point, let {Ak,J2f, Oj^) be the correspond- 
ing abelian variety with a (£n)-marking and denote by tt : A^ Bk the isogeny 
defined in Section[3l From the knowledge of {o-u)^^z(In)^ "-"^^ ''^'^ recover the co- 
ordinates the points of a maximal ^-torsion subgroup of B^ of rank g. Actually, 
even if ia.u)u^z(en) ^ valid theta null point it is possible to associate to 

('^f )tiGZ(£n}_ ^ .^-torsion points Pi G Bk- The main result of this section is 

Theorem |4] which states that a geometric point of Vj is non degenerate if and 
only if the corresponding Pi form a maximal subgroup of rank g of the -^-torsion 



points of Bk- To prove this Theorem, we introduce an action from the automor- 
phisms of the theta group to the modular space Mj;^. Using Theorem |4] and this 
action, we make exphcit the structure of Vj: we explain when two valid points 
give isomorphic varieties in Proposition [51 and how to obtain every degenerate 
points in the discussion following Proposition [71 

We start by making explicit the structure of the solutions of the algebraic sys- 
tem defined by J. For this let p : Z{n) x Z(l) — ^ Z{£n) be the group isomorphism 
given by {x,y) ^ £x + ny. Denote by Iq— the ideal of fc[y„|u G Z(ji)] for the 
theta structure On generated by the equations of Theorem [1] The homogeneous 
ideal Iq- defines a projective variety Vi^_, isomorphic to Bk- 

We have the following proposition [CL08j : 

Proposition 6. Let {o-v)^^z{In) ^ geometric point ofVj. For any i € Z(i) 
such that {ap(^j i))j^z{n) 7^ (Oi ■ ■ • jO), let Pi be the geometric point, o/P^^"'' with 
homogeneous coordinates io,p{j,i))jez{n)- Then for all i G Z(£) such that Pi is 
well defined, Pi is a £-torsion point ofViQ_. 

The proof of the preceding proposition in jCL08| proves moreover that if we 
denote by S the subset of Z(l) such that Pi is well defined for all i & S, then S 
is a subgroup of Z(I), the set {Pi, i G S*} is a subgroup of the group of ^-torsion 
points of V/q_ and the application i (z S Pi B[£] is a group morphism. 

Suppose that {0'v)j,^z(7n) ^ valid theta null point. Let {Ak,J^, Oj^) be the 
corresponding abelian variety with a (£n)-marking and denote by tt : Ak B^ 
the isogeny defined in Section [3l We can consider Ak as a closed subvariety of 
P^*-^"-* via the morphism provided by O-^. Using the action ([2]) of the theta 
group on iav)^^z(7n)^ '^^^ ^^^^ that for i G Z{£), the points with homogeneous 
coordinates {o,v+ni)y^z{ln) form the isotropic (for the commutator pairing) £- 
torsion subgroup Ki of Ak (with the notations of Section [3]). By definition of 
the isogeny tt, we have '^iiO'v+ni)yfzz(tP[)) — Pi and as a consequence tt{Ki) = 
{^Pi,i G Z{£)^. We see that if {av)^^z(hi) ^ valid theta null point then the 
{Pi)i^z(i) yfeli defined projective points which form a maximal subgroup of 
rank g of Bk[i]- Moreover, since the kernel of tt is K2, tt{Ki) — {Pi,i G Z(£)^ is 
the kernel of the dual isogeny tt : Bk ^ Ak- 

z(hi) is ^ general solution, it can happen that certain of the Pi are 
not well defined and as a consequence {o,v)y^z(7n) ^'^^ ^ valid theta null point. 
But even if every Pi are well defined, {civ)^^z{7ii) need not be a valid theta null 
point. We need a criterion to identify the solutions of J which correspond to valid 
theta null points. From the discussion of the preceding paragraph, we know that 
a necessary condition for a solution io-v)y^z{ln) '^^ >/ to be a valid theta null point 
is that (Pi)i^z(e) ^^'^ valid projective points which form a subgroup of rank g 
if Bk[£]. The Theorem [H asserts that this necessary condition is indeed sufficient. 
In order to prove this theorem, we have to study how a theta null point vary 
together with a change of the theta structure. 



We denote by Aut Ti.{S) the group of automorphisms ip of Ti.(S) inducing the 
identity on Gm,k- 



■ G^.fc ^ n{6) ^ K{S) 



^ G,„_fc ^ n{S) ^ K{S) ^ 

Obviously, the set of aU theta structures for ^ is a principal homogeneous space 
for the group Aut'H{d) via the right action 0s -ip = 0s o ip for -0 e Aut'H{S) 
and 0s a theta structure. So we can identify Aut7i((5) with the group of au- 
tomorphisms of theta structures. If ip is such an automorphism, it induces an 
automorphism ip of K{d). Denote by Sp{K{6)) the group of symplectic auto- 
morphisms of K{d). The preceding diagram shows that ip is symplectic with 
respect to the commutator pairing. Conversely, if ip G Sp{K{6)), we get an el- 
ement of Aut7i((5) defined by ip : {a,x,y) {a,ip{x),ip{y)). So the morphism 
•P' : Aut7i((5) Sp{K{6)), ip i-^ ip has a section that we denote by a. The kernel 
of If" is given by automorphisms preserving a symplectic basis and are determined 
by a choice of level subgroups Ki and K2 over maximal isotropic subspaces Ki 
and K2- It is well known [BL04, pp. 162] that such choices are in bijection with 
elements c € K{5): we map c G K{6) to the automorphism of 7i((5) given by 

{a,x,y)\-> {aes{c,x + y),x,y). (11) 

As a consequence, we get a split exact sequence 



^K{5) ^AutH(5) ^Sp{K{5)) ^0- (12) 

Suppose that 0s is symmetric, an automorphism ip g Aut7i((5) is said to be 
symmetric if it commutes with the symmetric action {a^x^y) 1— > (a, — x, —y) on 
^{{d). We denote by Autg H((5) the group of symmetric automorphisms of TL{5). 
Obviously, an automorphism ip S Aut7i((5) coming from c G K{5) is symmetric 
if and only if c g _fi'(5)[2] the subgroup of 2-torsion of K{5). 

Now consider (A^, 6*5) an abelian variety with a (5-marking and let {'&i)i^z{5) 
be the associated basis of global sections of .5f. Note that if ■(/; is a symplectic 
isomorphism of K{5) then ip = a{ip) is symmetric. We suppose that ip{Z{S))) = 
Z'l' X i'^, where Z'^ c Z{S) and Z'l' c Z{S). Denote by i^^)^ezi5) the basis of 
global sections of .if associated to {Ak,.^,0s-ip)- In the following, we give an 
explicit formula to obtain {'&i)i£Z{S) from the knowledge of {'di)i£Z{S)- 

Let ~ Ak/0s{ip{Z{S))) and vr : Afc — > A^ be the canonical map. The data 
of the maximal level subgroup 0s{iPii^T(^,y)y^z(S))) equivalent to the data 
of a line bundle ^0 on A^ and an isomorphism 7r*(^o) ^ Let Sq be the 
unique global section of Jffo, we can apply the isogeny theorem [Mum66[ Th. 4] 
to obtain 

#o = A7r*(io)- ^ (13) 



for Xek*. 

Now by definition we have 

=^((l,i,0))A. (14) 

wliere tlie dot product is the action ([2]). 

By evaluating at the basis of global sections of ^ in IT^ , we get an explicit 
description of the action of Sp{K{S)) on the geometric points of Ms- Actually 
the obtained formulas give a valid action of Sp{K{S)) on the geometric points 
of Ms- 

Now, let {o-u)u<£z(en) ^ geometric point of Vj. As Autg H{£n) acts on M-^, 
we are interested in the subgroup f) of Autg H{in) that leaves {0'u)uez(hl) ■ 

Lemma 1. Letip G Aut7i(£n). We say thatip is compatible withTi.(n) if it com- 
mutes with the morphisms p and i from ([6|) and ([7]) . Then 9) is the subgroup of 
compatible symmetric automorphisms ofT-L{tn). In particular it does not depend 
on {oLu)u£z(In) *^ '^^^^ subgroup o/Auts7Y(fe) that leaves Vj invariant. 



Proof- Let (A, Oj;^) be a triple corresponding to the theta null point io-u)u^z(h 
Let Tp e AvLtsTi{£n), and ia'u)uez{7^) = '^■("■u)uez(hl)- The Proposition [5] shows 
that {0''u)u^z(en) '^^^Y the associated theta structure O-^.ijj is 

compatible with the theta structure On of B. But this means exactly that ip is 
compatible with Ti.{n). 

We can describe more precisely the action of 9j: 
Proposition 7. The action of on Vj is generated by the actions given by 

where ip2 is an automorphism of Z{£n) fixing Z{n) and 

K)«GZ(i;i) ^ fee V'l -an) „g2(te)' (16) 

where ipi is a "symmetric" morphism Z{£n) Z{£) C Z{£n) and e-^ is the 
commutator pairing on TL{£n). 

Proof. Let ip G S). Since the exact sequence from equation (|12p splits, we only 
have to study the case where ip comes from a change of maximal level structure 
and the case where -0 comes from a symplectic base change of K{S). In the former 
case, let c € K{6) defining the symplectic base change by pTjl . Then c e ii'(J)[2] 
since ip is symmetric and from the compatibility conditions c G ■ip{Z(£)). As £ is 
odd, we have c = 0. 

In the latter case, ip can be represented in a basis {v^, VK)Ke{i g} of Z{£n) x 

Z(hi) by a matrix M[A, B,C,D]= (^^ ^ ^ e SP^^{Z). Since 

K ^Oj^i^izm cej^{z{hi)), 



we have i? = 0. SoD = *A ^ and we see that the action of is generated by 
the matrices 

1. M[A, B, C, D] such that C—0. Then A is an automorphism and the compati- 
bihty condition imphes that it must fix Z(n). Using ([TS]) and this yields 
the action 

2. M[^, B, C, L»] such that A = Id. Then *C = C. For x e we can write 
tp{{x,Q)) = (x, V'i(a;))- By looking at the conditions ^ and ([7]) we see that 

^{{x, y)) - (x, v) e ^{Z{1)) C (17) 

for all [x, y) E Z*(In) x Z*{In). Using dTT]), we deduce that ^1(2:) is in Z(l). 
In this case, we obtain the action (|16p following (fT3|) and p4|) . 

This completes the proof of the proposition. 

Remark 2. The action (jlSp gives an automorphism of the {Pi)i^z(I) "'^hile the 
action leaves the {Pi)i^z(e) invariant. In fact by taking a basis of Z{£n), we 
see that if ^ is a (£n)*''-root of unity, the actions are generated by 

where (aij")i,je[i,g] is a symmetric matrix and a^j £ li/nZ C Z/£n7j (via a; 1— > i?a;) 
for z, j e So each coefficient of one Pi is multiplied by the same ^'''-root of 

unity. 

From the preceding remark, we see that 9j leaves Vj invariant. Now, let 
^2 be a morphism of Z{£n) fixing Z(n). Here we do not require 'ip2 to be an 
isomorphism. We let -02 act on Vj by 

("")uez(£;i) ^ (^■>p2{u))uez(K) 

Since ip2 fixes Z{n) C Z{£n), it fixes the 2-torsion points in Z{£n), and it is easy 
to see that {o-ip2(u))u£Z(in) satisfies the equations of Theorem[2]and the symmetry 
relations. As a consequence, the point {0'4!2(u))uez(en) ^Jn- Moreover, as 
V'2 fixes Ziji), {(''■4,2{u))uez{In) ^ point in Vj, so we have a well defined action 
extending that of the form (|15p . 

By acting on Vj with a morphism of Z{£n) fixing Z{ri) which is not an 
isomorphism, we obtain a point of Vj which is degenerate: it is a theta null 
point such that the associated points Pi from Proposition [5] are well defined but 
not distinct projective points (so they do not form a rank g ^-torsion subgroup 
of Bk). 

There is another way to obtain degenerate theta null points in Vj . Take any 
geometric point {0'u)u^z(in) ^ ^'^'^ ^ subgroup 5 of Z(£) (in particular S is 
not empty). We define a new point {0''u)uez(ln) where 

^^■''^'^ I otherwise. 



Since i is odd, it is easily seen that {0''u)u£Z(7n) general a degenerate point 
in Vj: the Pi from Proposition [6] are not defined when i ^ S. 

Now, we explain that combining the two methods described above, we obtain 
all the degenerate theta null points of Vj. For this, let {o.'u)u£Z{hi) ^ degenerate 
point of Vj. Let S <Z Z(l) be the subgroup where the points of ^-torsion Pl,i & S 
of Proposition [6] are well defined. The points f / form a subgroup S' of the £- 
torsion points of Bk, and / : S' — > S",« i— > P/ is a group morphism (which 
may not be an isomorphism, since as io-'u)uez(ln) degenerate the P- are not 
necessarily distinct) . Now, we embed S' into a maximal subgroup T of rank g of 
and extend / to a morphism f : Z(i) ^ T (for instance if i ^ Z(i) then 
send i to the neutral point Pg). We take an isomorphism h between Z{i) and 
T. Theorem 2] that we prove later on shows that there exists a geometric point 
ifl'u)y^^z(In) ^ such that the corresponding group morphism i £ Z(i) i— > Pj 
is h. Now take ■02 to be the morphism of Z{£n) = Z(I) x Z(ji) which is the 
identity on Z(n) and h^^f on Z(£). Consider the point iaip2{u))uez(i^) with the 
coefficients p{j,i),i ^ S taken to be 0. Then it has exactly the same defined 
points P/ as Wu)uez{7n)- next lemma shows that it is the same point as 
i'^'u)u£Z(tn) ^'^ action of the form (fT6|) . 

We remark that the degenerate points in Vj are exactly the points where 
the action of ^ is not free: if iau)u£Z{m) ^ degenerate point such that the 
corresponding Pi arc not all well defined, then there is an action of the form 
giving the same point. If the Pi are well defined but do not form a maximal 
subgroup, then this time there is an action of the form ()15|1 giving the same 
point. 

By Remark [2] we know that if io-u)u^z{In) ^ theta null point giving the 

associated group {Pi, i G Z(£)}, then the points '4'-io,u)u£z(7n) where ip G Sj give 
the same associated group. In fact the converse is true: 

Lemma 2. Let icu)ii_^z{7n) '^'^'^ i^^)uez(In) ^^'^ geometric points ofVj giving 
the same associated group {Pi,i G Z(J)}. Then there exist ip E :^ such that 

(^")uez(i;r) — ■'/'•('^f )tiez(sr) ■ 

Proof. First, up to an action of type (jlSp . we can suppose that for all i G 
Z(i), we have p.''^"^"^^^*^^) ^ ^ rj.^^^ ^^^^^ ^^^^^ A, G fc such that 

{cp{j,i))jez{n) = A,(dp(j-,))je2(7,). Since ic^,)^^z(J^^) ^nd {du)^^z(I^) are projec- 
tive, we can assume that Aq = 1. We will show that up to an action of type 
for every i € Z{i) such that Pi is well defined, A^ = 1. But first we show that for 
such points, we have \[ — 1. 

Let i e Z(?) be such that {cpf^j,i))j(zz{n) is a well defined projective point. 
Let x,y,u,v e Z{2n) which are congruent modulo Z(ri), we remark that for 
/Li e {1, ...,£}, p(x, /u.i), p(y, i), p{u, 0), p{v, 0) are elements of Z{2£n) congruent 



modulo Z{£n). Calling Theorem [21 we obtain that 

( ^ X{t)Cp(x+y+t,{p+l).i)Cp{x-y+t,{tJ,-l).i))-{ ^ X{t)Cp{u+v+t.O)(^piu-v+t,0)) = 
teZ{2) teZ(2) 

— ( ^ X{'t)Cp(x+u+t,p.i)Cp(x-u+t,p.i))-{ ^ Xi't)Cp(y+v+t,i)Cp(y-v+t,i))j 
teZ{2) tGZ{2) 

(18) 

for any x G -^(S)- 

We have a similar formula involving {du)^^z{7^) - Using equation (fTS]) and an 
easy recurrence, we obtain that A^.i = A"^ where (up) is a sequence such that 
uq — 0, ui = 1 and it^+i + = 2.u^ + 2. The general term of this sequence 
is Up = ■ For ^ — £, we have 

Af = A,., = Ao = 1 (19) 

Now, by the symmetry relations, we have for j e Ziji), Cp(^j p i) = Cp(^_j^_p iy 

iplying this fc 
with (flQl) gives 



Applying this for n = 1 and j — 0, we obtain that Ai = A^ which together 



Af = 1 (20) 

which concludes the claim. 

Let (ei, . . . , eg) be the canonical basis of Z(I). Up to an action of type (fTB]) 
we may assume that Ae^ = 1 and Xa+ej — 1 for i,j G {1, . . . ,5}, j < i. Now let 
a,b G Z(I) be such that Aq = 1, Af, = 1 and Xa-b ~ 1- Then by Theorem[2]we 
have the relations: 

( ^ X{'t)(^pix+y+t,a+b)Cp{x-y+t,a-b))-{ ^ X{t)'^p{u+v+t,0)'^p(u-v+t.,0)) = 

tez(2) tez(2) 

= ( ^ X(i)Cp(a;+M+t,-6)Cp(a;-ti+t,6))-( ^ x(OCp(y+t)+t,a)Cp(y_«+t,a)) • 
teZ(2) teZ(2) 

(21) 

Since by symmetry, A_b = 1, the relations (PT|) give that Aa+b = 1. An easy 
recurrence shows that for any i € Z(I) we have A^ = 1, which concludes the 
proof. 

As a first application of this lemma, we have: 

Proposition 8. If i is prime to the characteristic of k and v 1 then Vj is a 
reduced scheme. 

Proof. We recall that Vj is the affine variety defined by J where J is the image 
of the homogeneous ideal / defining Mj;^, under the specialization map 



A;[a;ti|u S Z{£n)] k[xu\u G Z{£n),nu ^ 0], 



6ti , if u G Z{n) 



with {bu)u£Z{n) the theta nuU point associated to {Bk,^o,On)- 

By definition, Vj is a closed subvariety of the afhne space A'^^^"-'. For A G 
Z(£), denote by tta : A^^^"^ A^^") the projection deduced from the inclusion 
ipx : k[xu\u S Z{n)] k[xu\u S Z{in)], Xu ^ a;p(«,A)- In order to prove that Vj 
is a reduced scheme it is enough to prove that for any x geometric point of Vj 
and all A e Z(i), ii\{x) is a reduced point of A'^'"^. We consider two cases. 

If iT\{x) is not the point at origin of A^^"-' then it defines a projective point 
of P-^(^) which is a ^-torsion point of Viq_ by Proposition [B) As a consequence, 
'n\{x) is contained in the reduced line L between the origin point of A^^"^ and 
tta {x) . By the preceding lemma, the intersection of Vj with L is contained in a 
variety isomorphic to Spec(fc[x]/(x^ — 1)) where — 1 is a separated polynomial 
as £ is prime to the characteristic of k. We deduce that x is a reduced point of 
A^("). 

If TTx (x) is the origin point of A-^^") , it is enough to prove that tta (x) is reduced 
in the case that v — 2.ln fact, the set of equations generating J in the case v ^2 
contains the set of equations generating J in the case v = 2. We suppose now 
that V = 2. Let = (x„|u G Z{n)) be the ideal of k[xu\u G Z{n)] defining the 
reduced point at origin of A^("). Let Ja = J n ipx{k[xu\u £ Z{n)]) and denote 
by Ja«p the local ring of Ja in As J is a 0-dimensional ideal, we know that 
there exist m a positive integer such that Ja<p 3 in k[xu\u G Z{n)]r^. Let 
r\ be the smallest integer with this property. We want to show that ta = 1. In 
order to do so, we are going to use another formulation of the Riemann relations 
given by Theorem [1] 

For this, we let H(hi) = Z(hi) x Z{2) and H{n) ^ Z{n) x Z{2). We denote 
by p' : H{n) x Z(T) H{£n) the natural isomorphism deduced from p. For all 
V = {v',v") G H{in), we let yv = Y,tGZ(2) v"{t)xv'+t- Let ai, a2, as, 04, r G H{n) 
such that 2t ~ 01 — 02 — 03 — 04. Set ai = /o'(oi, 2A), a2 — p'(a2, 0), = p'{a^, 0), 
04 — p'io-i, 0) and ri — p'{t, A) so that we have 2ti = ai — a2 — as — 04. We 
write r = (t',t") and let H{2) = {x G ff(£n)|x is 2 -torsion modulo Z(2) x {0}}. 
By applying IMum67[ formula (C") p. 334], we have the following relation in J: 



yaiya2yoi3ya4 



1 

'2a 



E 

teH{2) 



(r"+t")(2i')y 



where t = {f ,t") G H{2). 

By definition, for i = 2,3,4, if we write Oi = (0^,0"), we have ^q,. = 
Stez(2) {t)bai.^f As by hypothesis {bu)u£Z(n) is valid theta null points, by 
applying [Mum67[ formulas (*) p. 339], we obtain that for any o^ = (0^,0") G 
H{n) there exists G 2Z{n) such that X)tez(2) i^)^a'.+i3'.+t 7^ 0. As a conse- 
quence, for any choice of oi, we can find 02, 03, 04, and t G H{n) such that 
2t = oi - 02 - 03 - 04 and for i = 2,3,4, j/q, ^ J2tez(2) ai'(*)^+t ^ 0- (We 
can take for instance oi = 02 = 03 = 04 so that oi — 02 — 03 — 04 G 2H(n) and 
then if necessary add to 02, 03, 04 elements of 2Z{n) in order to have 2/(ai,o) 7^ 0.) 



As an immediate consequence, we obtain that 7r2A(a;) is also the origin point of 

Let be the smallest integer such that ^ r\ and 4|r^. We remark that 
ip2\{k[xu\u e Z{n)]) — ^[yp'(t,,2A)b £ -^(")]- Let M be a degree monomial 
in the variables yp>(v,2\)- If necessary, by multiplying M by a suitable non null 
constant, we see that M is equal to a product M' of r'^/'i polynomials given 
by the right hand of (|22p . These polynomials have degree 4 and are sums of 
products of monomials of the form yp>(v.\) (using the symmetry relations). We 
deduce from this that M' g ^P''^ and as a consequence M' G J\. But this means 
that M G and as M can be any degree r^/4 monomial in the variables 
yp'{v,2\), we have proved that J2A<p 3 ^''>-/^. 

Let m be an integer such that 2™ A = A in Z(I). Using the previous result 
and an easy recurrence, we see that if > 1 then r\ = r2m\ < r\ which is a 
contradiction. 

As a second application of Lemma [21 we have: 

Theorem 4. Let {'^u)u^z(In) geometric point of Vj. For any i G Z(J), let 
Pi be the geometric point, if well defined, of P^*-"'' with homogeneous coordinates 
{0'p(j,i))jez(n)- Denote by S the subset of Z{£) such that Pi is a well defined 
projective point for all i £ S . If K — {Pi,i G S} is a maximal £ -torsion subgroup 
of Viq_ of rank g then {0'u)^^z(ln) ^^^^ defined theta null point. In other 

words there exists {Ak,J^,Oj^) an abelian variety together with a (in) -marking 
with associated theta null point io,u)u^z(In)' 

Proof. Let Ak be the quotient of Bk — V7©_ by K and let tt : Bk ^ Ak be 
the canonical isogeny. As if is a subgroup of there exists an isogeny 

TT : Ak Bk such that [£] —no n. Let ~ Tt*{^o). We are going to show 
that there exists a certain theta structure such that the theta null point 
associated to {Ak,^,0j^) is {au)uez(IFi)- 

Let K{^o) — Ki{^o) X K2{^o) be the decomposition into isotropic sub- 
spaces for the commutator pairing induced by the theta structure 0^. Denote 
by K the kernel of tt. As ^ — 7r*(^o), we know that K is an isotropic subgroup 
of K{^) for the commutator pairing. Moreover, by construction it is contained 
in the ^-torsion subgroup of A and by hypothesis has rank g. We choose a de- 
composition as isotropic subspaces K{Ji') = Ki{J^') x K2{^) such that K is 
contained in if2(if) and for i = 1,2, Tr{K,{^)[n\) = K,{^o). 

Denote by k : G(^) K{^) the natural projection. By the descent theory 
of Grothendieck, there exists a unique level subgroup Ki of G{^) contained in 
K~^{K) such that the quotient of (Ak,^) by the action defined by Ke gives 
(Bfc,^o). Let G*{^) be the centralizer of Ki in G(^). By |Mum66l Prop. 2 pp. 
291], we have an isomorphism 



Let G(if)[n] = K-^{K{.^)[n]). We remark that 



1. G(^)N is contained in G*(^), 

2. n{G{J^)[n] n K() is the zero subgroup of Ak- 

Let K^) be the level subgroup of G(^o) defined as the image by On of the 
subgroup {^T^Ty)y,^z{n) of'H{n). An immediate consequence of 1. and 2. is that 
there exists a unique level subgroup Kn of G(.5f) such that i{Kn) = -K'q. 

Denote by K2 the level subgroup of G(^) whose restriction over K{^)[n] 
and if (^)[^] is respectively given by Kn and i^T^. By construction, we have 

i{k2) = ko. (23) 

Choose any theta structure 0-^^ : Ti{ln) G(J£) such that the image by 
of the subgroup {X,^-,y)yizz(tn) exactly if2- Because of ([23l) and construction 

of Proposition [H we have zJq'^" = 7r*(i?^"). 

We suppose moreover that Oj^ is such that for all x E Z{n), i{0j^{l, x,0)) = 
6*71(1, X, 0), where we consider Z{n) as a subgroup of Z{in) via the map x t-^ £x. 
We remark that by construction, Oj^ and On verify the conditions (O and ([7|) 
and as a consequence are 7r-compatible. As a consequence of Corollary [1] we 
have that for aU ie Z{n), df^ = TT*{^f^). 

Let {a'u)uez(in) ^'^ theta null point associated to {Ak,^,Oj;^). For i G 
Z(n), denote by Qi the geometric point of P^'"-* with homogeneous coordinates 
("p(j,i))je2(")- 

We know that the projective coordinates of a maximal isotropic i?-torsion 
subgroup of Ak is obtained by the action of the theta group on io-u)u£Z(i^) 
translation. Denote by K' the ^-torsion subgroup of Ak given by the points with 
projective coordinates {0'v+i)vez(ln)- construction, K' is the dual of K for 
the commutator pairing which implies that Ak is exactly the quotient of Bk by 
7t{K'). As a consequence, we have tt{K') ~ K. 

The applications Z{1) — > Bk[P\, j 1— *■ Pj is a group morphism (see for instance 
the proof of [CL08[ Lemma 5.6]), as well as the application Z(l) Bk[i], j n- 
niQj). By changing the theta structure Oj;;^, we can suppose that for all j e Z(I), 
TT{Qj) = Pj. As a consequence, for j S Z(I) there exists Xj S k such that for 
iGZ{n), ap(j^i) = Aja^(j. We know moreover that {au)^^z{hl) and «)„gz(te) 
are geometric points of Vj. Applying Lemma [2] we are done. 

If {o,u)uez(I7i) is a geometric point of Vj, we denote by G{.{.aii)u(iz(hl)) t^'^ 
subgroup of Bk[P\ generated by the valid projective points {ap(^ji))j^z(n) for 
i E Z(l) of Viq_ = Bk- The preceding theorem tells us that whenever a solution 
{o-u)uez(tn) of is such that G{{au)^^z(In)) ^ maximal i!-torsion subgroup 
of Bk[£] then it is a valid theta null point, that is, it corresponds to a certain 
(Ak , 1 Oj^) . It would be desirable to be able to determine which maximal rank 
g subgroups of Bk[£] can arise as a G{x) where a; is a geometric point of Vj 
representing a valid theta null point. 

For this, let ^0 = [^]*-2o on Bk- As J^'o is symmetric, we have that ^0 — 
Ji'g and as a consequence K{^o), the kernel of ^0 is isomorphic to Z{Pn)- 



The polarisation induces a commutator pairing ej(^ on K{^q) and as 
descend to via the isogeny [£], we know that ej(„ is trivial on Bk\P\- For 
x\,X2 G Bk\(\, let x'ljx'a € Bu\f'\ be such that f.x- = Xi for i = 1,2. We 
remark that x\ and defined up to an element of \i\ ■ As a consequence, 

e^o(a^i,a;2) = e_^p(a;i,X2), does not depend on the choice of x\ and x'2 and if 
we put ey/{x\TX-2) = (x']^, 0:2), we obtain a well defined bilinear application 
ew '■ Bk[i] X Bk[l] — > fc. As e is a perfect pairing, for any x[ G Bk[£^] there 
exits 2:2 G Bi^if^] such that e ^^(x'j^jXj) is a primitive root of unity. As a 
consequence, for any xi G Bk[l] there exists 2:2 G B]^[£] such that en/(a;i,a;2) is 
a primitive root of unity and ew is also a perfect pairing. 

Theorem 5. A maximal £-torsion subgroup of Bk of rank g is of the form G{x) 
where x is a geometric point of Vj corresponding to a valid theta null point if 
and only if G{x) is an isotropic subgroup for the pairing ew- 

Proof. Let {au)u£z(in) ^ geometric point of Vj corresponding to a valid 
theta null point. We know that {(iu)u^z(7n) theta null point of a triple 

{Ak , -Sf, O-^) . The theta structure 0-^ induces a decomposition -ft'(^) = Ki {^) x 
i^2(^) into isotropic subgroups for the commutator pairing e_sf . As the isogeny 
TT is such that 7r*(z9f") = i^f^" for all i G Z{n) (and identifying i e Z{n) 
with £i G Z{£n)), we know that G{{au)^^z(M)) — 7r(-fCi(^)). We denote by 
TT : Bk ^ Ak the isogeny such that tt o tt — [£] a,s in the diagram For any 
Xi,X2 G G{(au)^^z^j;^'f), there exists a;i,a;2 G Ki{^)[£] such that Xi = 7r(a;i), 
i — 1,2. Let x'j e be such that £.x'i = xi. We have ew{xi,X2) = 

e^o{x'i,X2) = e^(7r(a;i),7r(a;2))- But 7r{x2) = n o tt{x2) = [^1(^2) = 0. As a 
consequence, we have ewi^i, X2) = 0. 

Now, we prove the opposite direction. Let G be a maximal rank g ^-torsion 
subgroup of Bk [£] which is isotropic for the pairing ew and G be the dual group 
of G for the pairing ew- As ew is a perfect pairing, G is also a maximal rank g 
^-torsion subgroup of Bk[£]- We want to show that G is of the form G{x) with x 
a geometric point of Vj where J is defined by the triple {Bk,.^o,0n)- For this, 
we consider the isogeny n : Bk A^ with kernel the subgroup G oi Bk- As G is 
contained in Bk[£], G is an isotropic subgroup of {Bk,^o), and ^0 descend via 
TT to a polarization ^ on Ak- Let tt : Ak — > Bk be the isogeny with kernel tt{G). 
By the commutativity of the following diagram, 




^ descends via tt to .2o. 



The theta structure 0^ induces a decomposition K{^i^) — A'i(^o) x K2{^^)- 
Let Xi — Tr{x'j) with a;^ G G and « = 1, 2. Let y'^ E Bk[£'^] be such that £.y[ — x'^. 
We have by hypothesis 1 = ewix'i,x'2) — e^o(2/'ii^2) ^^'^ a consequence 
1 = 6^0(0;']^, = e^{xi,X2)- Thus 7r(G) is isotropic for the pairing e^. As 
a consequence, we can chose a decomposition i^(-Sf) = Ki{^) x iir2(-S^') such 
that for i = 1,2, tt{K,{^)) = K,{^o) and i^2(^)[^] = 7r(G). Take any theta 
structure 0-^ for ^ compatible with this decomposition. Let {o-u)u^z(en) 
associated theta nuU point. By Corollary[Tl io,ii)j^^z{7n) ^ geometric point of 
Vj. Moreover, we have G{{au)^f=z{Wi)) = = G. 

Our study of vaUd theta null points allows us to better understand the 
geometry of Vj. We know from Proposition O that Vj classifies the isogenies 
TT : Ak Bk between marked abelian varieties verifying the compatibility con- 
dition. 

Taking the dual of tt gives an isogeny from Bk to Ak with kernel K = 
Tr{Ki) — {Pi, . . . , Pji}- Thus the theta null points on Vj correspond to varieties 
Z-isogeneous to Bk- But we have seen in Proposition [7] that it may happen that 
different points of Vj give the same kernel K and hence the same isogeneous va- 
riety. We want to classify the points of Vj corresponding to isomorphic varieties 
£-isogeneous to Bk- 

To do that, let K he a maximal isotropic subgroup of rank g of the points of 
^-torsion of Bk- We are interested in the class 1^ of isogenies of kernel K. More 
precisely, if tt is an isogeny from Bk to Ak with kernel K, then Ik is the class 
of isogenies tt' : Bk ^ A'f. such that there exists an isomorphism t/j : Ak A'f. 
that makes the following diagram commutative: 




Proposition 9. Let K be a maximal subgroup of rank g of the points of (.-torsion 
of Bk which is isotropic for the pairing ew- There is a point {0'u)ui^z(ln) ^ 
such that the corresponding dual isogeny tt : Bk — ^ Ak is in Ik- Every other 
point in Vj giving the class Ik is obtained by the action of Sj on {o^u)u£Z(In)- 
In particular, the geometric points ofVj/Sj are in bijection with the -isogenies 
ofB- 

Proof. Let K = {Pi,i E Z(l)} be such a maximal subgroup. Theorem [4] gives 
a geometric point {0'u)y^^z(en) corresponding to a marked abelian variety 

{Ak,.^A,0A) such that the associated isogeny n : Ak ^ Bk sends Ki{^a) to 
K. Hence, the unique isogeny n : Bk ^ Ak such that n o n = [^], is in Ik- If 
i'^'u)uez(en) another valid theta null point in Vj, corresponding to a marked 



abelian variety [^A! ^^a' ^ ©A') such that the dual of the associated isogeny gives 
the same class as tt, then we have the following diagram: 




By definition of the associated isogenics tt and tt', we know that ££a = tt*{^b) 
and = tt'*{^b) = iP*{-^a)- So tp induces a morphism of the theta groups 
G{^a) and G{^a'), and puUing back by the theta structures we get a symmetric 
automorphism ■0 of Ti{in). Since the theta structures Oa and 0a' are compatible 
with 6>B, %() is in ij. This shows that (au)„g2(£;^) and {0''u)iiez(ln) ^^'^ same 
orbit under ^. 

Together with the study of degenerate theta null points, it is now possible 
to count the points in Vj. For instance, take 5 = 1, n = 4 and i — i. Let E 
be an elliptic curve, and (&u)ti6(z/riZ) be a level 4 theta null point on E. There 
are 4 = #P^(F3) classes of 3-isogenies from E^ and 6 = 3 x 1^9(3) solutions in 
Vj for each class. The actions (US]) are given by {au)u<^(i./tni) ^ {ax.u)u&z/in-L 
where x € TLjtnL is invertible and congruent to 1 mod n. There are <p(i') such 
actions. The actions are given by (au)u&linL ^ iC'" a«)uez/£nZ where C 
is a ^*''-root of unity and c G Z/£Z. 

If g = 2, it is easy to compute the number of valid theta null point in Vj. First, 
we remark that the number of isogeny classes of degree i'^ of a given dimension 2 
abelian variety is parametrised by the points of a Grassmanian Gr{2,4:){¥g) 
which are isotropic (see Theorem[5]): there are {£'^ + 1){£ + 1) such points. 

Next, the number of actions of the form (jlSp is parametrised by the number 
of invertible matrices of dimension 2 with coefficients in with is given by 
(£^ — 1).(£^ — i). The number of actions of the form (fT6|) is £^ (the number of 
symmetric matrices of dimension 2). As a consequence, the number of valid theta 
null point in Vj is 

We remark that this number is a 0{£^^). For g ~ 2, £ = 3, we have 51840 valid 
theta null points in Vj. 

For a general g and £, we assess the order of the number of valid theta null 
point which arc solution of Vj. The number of isotropic points of a Grassmanian 
Gr{g,2.g)i¥i) is a 0(^9(s+i)/2). The number of action of the form HSl) is a 
0{£^ ) and the number of action of the form is a 

0(£s(s+i)/2). We deduce that 
the number of valid theta null point in Vj is bounded by 



0(£2-9'+9). 



(25) 



Example 1. In the case of genus 1 and small £ it is possible to list all the solutions 
of Vj. We take t and let E be the elliptic curve given by an affine equation 
2/2 = x'^ + ll.a; + 47 over Fyg. A corresponding theta null point of level 4 for E 
is(l:l:12:l). The four subgroups of 3-torsion of E are 





= {(1 


1 


12 


1),(37 


54 


: 46 


:1),(8: 


60 : 


74 : 1)} 


K2 


= {(1 


1 


12 


1),(67 


10 


: 68 


:1),(62 


: 8 : 


70 : 1)} 




= {(1 


1 


12 


1),(42 


5 : 


15 : 


1),(40: 


16 : 


3:1)} 


Ki 


= {(1 


1 


12 


1),(72 


56 


: 31 


: 1),(69 


: 24 


: 33 : 1)} 



All geometric points of Vj are defined over F79(u) where w is a root of the 
irreducible polynomial + 9.X + 76. For each of the four subgroups Ki, there 
are 6 geometric points of Vj giving the curve E/Ki. We give a point in each 
class (the other points can be obtained via the actions (fT5|) and P^ ): 

Qi = (16u2 + 19u + 17 : 1 : 46 : IBv^ + 19u + 17 : 37 : 54 : 34z;2 + 70u + 46 : 
54 : 37 : 16^2 + 19w + 17 : 46 : 1) corresponds to Ki. 

Q2 = (64^2 + 67u + 68 : 1 : 68 : 64^2 + 67u + 68 : 67 : 10 : 57^2 + 14u + 26 : 
10 : 67 : 64^^ + 67u + 68 : 68 : 1) corresponds to K2. 

gg = (8i;2 + 49w + 48 : 1 : 3 : 8^2 + 49w + 48 : 40 : 16 : 17^^ + 35w + 23 : 16 : 
40 : 8v^ + 49w + 48 : 3 : 1) corresponds to K3. 

Qi = (32^2 + 73u + 34 : 1 : 33 : 32^2 + 73w + 34 : 69 : 24 : 681;^ + 7u + 13 : 
24 : 69 : 32^2 + 73u + 34 : 33 : 1) corresponds to Ka. 

We also have the following degenerate points in Vj: if we take x = 9 in the 
action (fT5)) . the image of the class of any is C = {(55 : 1 : 12 : 55 : 1 : 1 : 28 : 
1 : 1 : 55 : 12 : 1), (1 : 1 : 12 : 1 : 1 : 1 : 12 : 1 : 1 : 1 : 12 : 1), (23 : 1 : 12 : 23 : 1 : 
1:39:1:1:23:12:1)}. For this class, the corresponding ^-torsion subgroup 
(the points P, of Proposition O is {(1 : 1 : 12 : 1), (1 : 1 : 12 : 1), (1 : 1 : 12 : 1)} 
which has rank 0. On C the action is trivial, so there are only 3 points in 
this degenerate class, coming from the action (|16l) . The last degenerate point is 
(1 : : : 1 : : : 12 : : : 1 : : 0), alone in its class. 

We conclude this section with some remarks concerning the case v = 1 and 
the case where the characteristic of k is equal to I. First, for computational 
reasons, for instance in order to limit the number of variables when computing 
the points of Vj, we would like to have v as small as possible. All the results of 
Section [5] are valid under the hypothesis that v ^ 2 and that the characteristic of 
k is different from i. In the case = 1, we can not even prove that Vj is a zero 
dimensional variety. Nonetheless we have made extensive computations which 
back the idea that even in the case = 1, in general, Vj is a zero dimensional 
variety whose degree is of the same order with respect to the parameter H. as in 
the case v = 2. 

In the case that the characteristic of the base field k is equal to i and ^ 2, 
the proof that Vj is a 0-dimensional scheme is still valid. In this case Vj is not 
anymore reduced and the computation of the number of solutions of Vj are not 
valid. Nonetheless, from our computations, we see that in this case the degree 



of the variety Vj is of the same order with respect to the parameter £ as in the 
case where the characteristic of k is different from £. 

In the following section, we give an algorithm to find the solutions of Vj. We 
can prove that this algorithm is efficient in the case v ^ 2 and when the charac- 
teristic of k is different from £. In the case that v — 1 or when the characteristic 
of k is equal to £ we will make the hypothesis that Vj is a zero dimensional 
variety whose degree is given by formula ([^5]) . Under these hypothesis, we can 
also prove that our algorithm is efficient. 

6 An efficient algorithm 

We would like to use the formulas of Section 2] to compute the image of the 
modular correspondence for some positive integer £. We have seen that the 
main algorithmic difficulty is to solve the polynomial system defined by the 
equations of Theorem [T] together with the symmetry relations. The aim of this 
section is to give an algorithm to solve efficiently this system. We have made an 
implementation of our algorithm and used it to test the heuristics described at 
the end of Section [S] 

Let n = 2'^ . In this section, fc is a finite field. We let {Bk, ^o, On) be a dimen- 
sion g abelian variety together with a ri-marking and we denote by {bu)u£Z(n) 
its associated theta null point. Let J be the image of the homogeneous ideal 
defining Mj^ given by the equation of Theorem (TJ under the specialization map 



We denote by Vj the 0-dimensional affine variety (heuristically 0-dimensional 
if J/ = 1) defined by the ideal J. Let p : Z{n) x Z(£) Z{£n) be the group 
isomorphism given by {x,y) i-^ £x + ny 

6.1 Motivation 

In order to find the points of the variety Vj a first idea is to use an efficient 
Grobner basis computation algorithm |BW93j such as F4 |Fau99| . We have car- 
ried out computations in the case 5 = 2, 1/ = 1 and £ = 3 with respect to a total 
degree order (the DRL [AL94ICL092] or grevlex order) using the computer alge- 
bra system Magma |BCP97j implementation of F4. From our computation, we 
could conclude that 

— even for a small coefficient field (fc — F310), it takes 20 hours of computations 
using Magma on a powerful computer with 16 Go of RAM; 

— as expected from the computations of Section O the number of solutions 
in the algebraic closure fc of fc is big: 30853 solutions in characteristic 3 
(We note that this is coherent with the number of solutions discussed after 
Proposition [9] when g — 2, v — 2 and 1? = 3). 



A;[a;„|u G Z{£n)] k[xu\u E Z{£n), nu 7^ 0], 




— to fully solve the system (that is to say, find explicitly all the solutions in k) 
we need to compute a second Grobner basis with respect to a lexicographical 
order. 

This last operation can be done using the FGLM [FGLM93] algorithm. In our 
case it is equivalent to compute the characteristic polynomial of a 30853 x 30853 
matrix. This computation did not finish using Magma for the base field k — F310 . 
So we see that even for g = 2, 1/ = 1 and £ = 3 the computation of the points of 
Vj is painful using a generic algorithm. In this section, we give an algorithm to 
solve efficiently the algebraic system defined by J for small i over a big coefficient 
field. As an application of our method, we can mention the initialisation phase 
of a point counting algorithm |CL08| . 

The main idea of our algorithm is to use explicitly the symmetry inside the 
problem deduced from the action of the theta group: we compute a Grobner basis 
not for the whole ideal J but rather a Grobner basis of a well chosen projection 
Jn k[Xp(^y^x)\v € Z(n)] for A G Z(l). With our strategy, the same problem 
{k = F310) can be solved in seconds and far bigger problems {k = F31500) can be 
solved in less than 1 hour (see Section [6^ for experimental results). 

6.2 Assumptions 

Our method is a combination of existing algorithms. We first describe in full 
generality the assumptions upon which our algorithm is faster than a general 
purpose Grobner basis algorithm. Then, using the results of Section [5l we ex- 
plain that these assumptions hold for J in the case that > 2 and that the 
characteristic of k is not ^. If = 1 or if the characteristic of k is equal to £, we 
can not prove the assumptions but we have made extensive computations which 
show that in general our algorithm is much more efficient than a general purpose 
Grobner basis algorithm. 

Let T be a set [a;i, . . . , Xs] of variables, we assume that J C A:[r] is a zero di- 
mensional ideal generated by the polynomials [/i, . . . , fm] where for i = 1, . . . , m, 
fi is a polynomial in k[T]. We make the hypothesis that we can split the set of 
variables into two subsets T — X UY such that the ideal K = J O k[Y] contains 
low degree polynomials. 

In order to make precise what we mean by low degree polynomials, we denote 
by Igen an ideal generated by the polynomials [gi, . . . , gm] where for i = 1, . . . , m, 
gi is a general polynomial of total degree deg(/i). We define for any ideal / of 
k[T]: 

Dy{I) - min{deg{g) | ^ .g e ^ H k[Y]}. 
Our assumption that Jn k\Y] contains low degree polynomials means that 

Dy{J) « Z?y(/ge„). (Hl-1) 

The previous assumption implies that our algorithm will perform much faster 
with the particular ideal J than it would do for a general ideal Igen ■ 



We must also ensure that it is more efficient to compute a Grobner basis for 
J n k[Y] instead of a Grobner basis for J. If we suppose that a Grobner basis 
computation for a total degree order has the same complexity for J and Igem 
we have to check that Dy{J) ^ DT{Igen)- It is well known that, generically, 
a lower bound for DT{Igen) is given by the Macaulay bound which is given by 
DT{Igen) = 1 + X)i=i('i'2S(/j) — 1) if ™ ^ s- We can now state explicitly the 
second part of our first assumption: 

m 

i?y(J)«^dcg(/,). (Hl-2) 

1=1 

Our second assumption is that J can be decomposed into many prime ideals. 
There exists a positive integer r ^ 1 such that 

\/j ~ Pi n ■ ■ ■ n Pr and Pi is a prime ideal. (H2) 

We recall that for a homogeneous ideal we define the Hilbert function HF/ (d) = 
dim(k[T]/I)d and the degree of the ideal /, deg(/), is given by the Hilbert series 

^ HF/(d) z' = f.^'i^r) and deg(/) = Af(l) ^ 0. With this, we can state the 
1=0 ^ ^' 

third (optional) assumption 

deg (VT) < deg (/) . (H3) 

We discuss the validity of hypothesis (|Hl-ip . (|H1-2|1 . p!2|l and (|H3| in the 
case that J is defined as in the introduction of the present section. First, we 
remark that Dyilgen) can be easily computed: let M{s,d) be the number of 
monomials of degree less or equal to d in s variables. The total number of 
solutions counted with multiplicities of /gen is given by the Bezout bound: 
D = n^-^ deg(/i). Hence, we have 

DY{Igen) = min{M(/i, d) > D} , (26) 

d 

where h is the cardinal of Y and M{h,d) ~ {^^'^)- By considering M{h,d) as 
a polynomial in the unknown d, we obtain that for a given h, Dyilgen) is the 
biggest real root of the polynomial: 

■ 1=1 

As a consequence, we have 

Dyilgen) -D^ooih\D)i. (27) 

We know moreover that Mj:^ has dimension 1/2. (5 +1) and is embedded via 
the relations given in Theorem [T] in the projective space of dimension in£)^ — 1. 
We deduce that J contains at least (ni!)^ — 1/2.(7.(5 + 1) algebraically independent 



polynomials. As the equations of Theorem [T] have degree 4, a lower bound for D 

is /^(niy-l/2.g.(g+l)-l^ 

On the other side, if we chose for j g Z{T), Y = [xp(^u,j)\u G ■^("■)]; "^"^ know 
by Proposition [5] that the solutions of the system J n k[Y] can be either the 
origin point of A'^^"^ or represent a ^-torsion point of Vig_. In this last case, 
by Lemma [2] we know that there is £ solutions of J corresponding to the same 
projective points. Denote by D' the number of solutions of J n k[Y] counted 
with multiplicities. We have D' ^ ^2g+i _|_ ^ ^^d using the heuristic evaluation 
of Dy{J) given by (pS)) . we obtain 

DY{J)-^D^oo{h\D')i. (28) 

For a fixed g and the cardinal of Y are fixed. Using (P7)l and ([^5]) . we see that 
hypothesis (|H1-1[) is verified for £ big enough. 

Next, 1 + E™i(deg(/i) - 1) = 3.(n€)9. On the other side, {h\D')T: with 
D' ^ £2g+i _|_ g^jjj /i = n^. As n > 2, we have, using the Stirling approximation 
formula, that {h\D')h = 0{£) and hypothesis (|Hl-2p is verified as soon as ,g ^ 2 
and £ big enough. 

Since we want to find at least one solution of J defined over fc, we can assume 
that such a solution exists. By Proposition [SI this implies that there exists a 
subgroup G of rank at least 1 of the ^-torsion group of Vi^_ such that all the 
points of G are defined over k. As the solutions of Jn k\Y] are points of V7q_[£], 
we conclude that for r ^ ^ we have: 

\/j = Pi n ■ ■ ■ D Pr and Pi is a prime ideal, 

and hypothesis (jH2p is verified. 

In general, we know from Proposition [8] that the hypothesis (|H3[) is not 
verified since J is a reduced ideal. Nonetheless, in the case that the characteristic 
of k is equal to £, the scheme defined by J is not reduced and we use l|II3|) in 
order to speed up the computations. 

6.3 General strategy 

In the following, we give a general strategy for computing the solutions of the 
algebraic system defined by J. All the steps of our algorithm are standard with 
the exception of step 1 and step 4. In step 1, we try to use as much as possible 
the assumptions (jHl-ll) and (|Hl-2p and step 4 is based upon the assumptions 
l|II2]) .(lH3]). 

Step 1 Using a specific algorithm given in Section [6. 4| we compute a truncated 
Grobner basis for an elimination order and a modified graduation. This 
allows us to obtain an zero dimensional ideal Ji contained in J. In gen- 
eral Ji is not equal to J. The output of the algorithm is a sequence of 
polynomials [pi, . . . in k [Y] such that Ji is generated by (pi, . . . ,Pk)- 

Step 2 Compute a Grobner basis Gdrl of Ji for a total degree order (DRL or 
grevlex). This can be done with any efficient algorithm for computing 
Grobner basis, for instance F4. 



Step 3 Compute a Grobner basis Glcx of Ji for a lexicographical order. This can 
be done by using the FGLM algorithm to change the monomial order of 
Gdrl • 

Step 4 Compute a decomposition into primes of the following ideal: 

yjj" = Pi n • • • n 

We assume that deg(Pi) ~ 1 (if it is not the case we replace k by some 
algebraic extension of k). 

Step 5 For i from 1 to r, we repeat the following Steps a,b,c for the ideal {Pi)+I: 

(a) Compute a Grobner basis Gi of {Pi) + I for a total degree order 
(DRL). 

(b) Change the monomial order to obtain G[ a lexicographical Grobner 
basis of {Pi) + I . 

(c) Compute a decomposition into primes: ^fpA^ — n • • • n P,; 
(by convention j_i = r). 

bmce we have sjl = ^fT[T^I = y/P[ n / n • • • n ^/ Pr r\ I and since the decom- 
position of each component ^/Pi H / is done by step 5 of the previous algorithm, 
we obtain a decomposition of the ideal /: 

\/7 = P,+i r\---C]Pj^ 

Remark 3. Once we have obtained a point P of Vj corresponding to a valid 
theta null point, we can recover all the solutions of Vj corresponding to valid 
theta null points using the action given by Proposition [71 

6.4 Description of the algorithm 

In this section, we give a detailed explanation of the Step 1 and Step 4 of the 
algorithm described in Section [^751 
Step 1: elimination algorithm 

The normal strategy for computing Grobner bases (Buchberger, P4, P5) con- 
sists in considering first the pairs with the minimal total degree among the list 
of critical pairs (see |CL092IBec93| . for instance). 

In the following, to select critical pairs, we consider only the total degree 
with respect to the first set of variables X. More precisely: 

Definition 1. Partial degree 0/ critical pair p — {f,g): 

degx {p) — total degree of 1cm ^LT(/), LT(g)^ 

in the polynomial ring R [X] where R = k [Y] . 

Moreover, we stop the computation of the Grobner basis as soon as we find 
a zero dimensional system in k [Y] . Consequently we obtain an new version of 
the F4 algorithm: 



Algorithm 6 Algorithm F4 (modified version) 



!F a finite subset of k[xi, . . . , Xs] 
< a monomial admissible order 
X ^ [xi,...,Xf,] and Y = [x^+i, . .. ,Xs] 
Output: a finite subset of k[xi, . . . , Xg] . 
G:=F and P := {CritPair(/, g) | {f,g) £ with f ^ g] 
while P ^ and dim(G n fc [F]) > do 

d := min{deg^ (p) I P G P} minimal partial degree of critical pairs 
extract from P, Pd the list of critical pairs of degree d 
R :==MATRix_REDUCTiON(Left(Pd) U Right(Pd), G) 
for h G R do 

P-~PU {CritPair(/i, g) \ g G G} 
G -.^Gu lh} 
return G 



Step 4: decomposition into primes 

The known general purpose algorithms to compute a primary decomposition 
of an ideal are inefficient in our case. To speed up the computation, we proceed 
following the three steps: 

Step 1 The basis Glox always contains a univariate polynomial g{xs). We can fac- 
torize this polynomial. We will see that this is the most consuming part of 
the whole algorithm. We obtain 

g{x,) = f,{xsr----fi{xsr'. 

Step 2 For all factors i from 1 to Z we apply the lextriangular algorithm [Laz92j to 
obtain efficiently a decomposition into triangular sets of Ji + {fi (xs))- 
We can describe the algorithm beginning by the special case of two variables 
[xs-ijXg] (this enough in our case since we assume that fc = fc as we will 
see later). By a theorem of Lazard Laz85j Theorem 1], the general shape of 
Glcx the lexicographical order Grobner basis is as follows: 

hi {Xs^l,Xs) 
, h2iXs-l,Xs) 

hg (^s — 1 1 •^s) ~ -^s—l * * * 

^ polynomials in variables Xi, . . . ,Xs 

with ki < k2 < ■ ■ ■ < ks and gi (xg) \ g2 {xs) \ ■ ■ ■ ■ Hence we can obtain for 
free some factors of g{xs)- 



91 i^s) f <ii + 

92 (xs) tx';ti + 



(29) 



step 3 



9KXs) = , . gi{xs) 

For any factor fi {xg) of 5(2;^) = fi{xs)°'^ ■ ■ ■ fi{xs)°'', it is enough to find 
the first element hj Xs) of the Grobner basis such that 

gcd(/i ixs),gj (xs))^ 0. 

In our case k = k and each factor is hnear fi (xg) = Xg — f3i so that we search 
for the first j such that gj (Pi) 7^ 0: then we obtain a new polynomial in 
one variable hj {xs~i,Pi) that can be factorized. Hence we can iterate the 
algorithm for all the other variables Xs-2, ■ ■ ■ ,xi. 

6.5 First experiments and optimizations 

In this section, we give running times for an implementation of the strategy that 
we have presented in Section [6.2l We also explain some important optimizations. 

The main motivation of the examples presented is this section, is to show 
that the initialisation phase of the point counting algorithm described in [CL08| 
can be made efficient enough to be negligible in the overall running time of the 
algorithm. For this, we take g = 2 and n = 2 and we work over a field k of 
characteristic 3 or 5. We construct a theta null point of level 2 corresponding to 
an abelian variety Ak of dimension 2. We construct the modular correspondence 
of level £ where £ is the characteristic of k. Any valid solution of the modular 
correspondence will corresponds to the theta null point of level 2£ of an abelian 
variety isogeneous to Ak- We can then use the algorithm of |CL08| to count the 
number of points of Ak- 

First experiments As explained in 16.11 if we can try to compute directly a 
Grobner basis of the ideal generated by the equations, even when k is very 
small (fc = F310 for instance), it takes 10 hours of computations on a powerful 
computer with 16 Go of RAM just to compute a DRL Grobner basis. Moreover, 
in characteristic 3, there is a huge number of solutions: 30853. This imply that 
there is no hope to solve efficiently the corresponding problem directly. 

Keeping the notations of the beginning of Section [6l we apply the method 
described in 16.31 to find the solutions of J. We let = 1, £ = 3 and g = 2 
so that Z(hi) = (Z/6Z)2. Let T = [xu\u G Z(£^)]- For j € Z(l), we define 
Y = [a;p(„j-) |tt e Z(fi)]- Taking j ~ p(0, 1) and in the following, for u = G 
Z{£n), we let Xu = Xij. With these notations, we take Y = [0:31, ^32, X02, xqi] and 
X — T — Y the set of all other variables. Then we consider J embedded in the 



polynomial ring k[T] where k is F3fc or ¥^k . In that case JOk [x^i, 2:32, X02, a^oi] = 
J n fc [y] is an ideal of degree 160 (to be compared with 30853 the degree of the 
whole ideal J). When k — ¥^k (resp. k — W^^k) the polynomial g(xs) obtained 
in section [6.41 is a square- free polynomial of degree 124 (resp. a non square- free 
polynomial of degree 70). We report in the following table some first experiments 
using the algorithm of section 16.31 implemented in Magma and in C (see section 
16.61 for a full description of the experimental framework). First we consider only 
very small example: 



AlgolOl l Step 1 Step 2 Step 3 Step 4 Step 5 

k = F510 0.35 sec 0.25 sec 3.24-1-0.01=3.25 sec 8.0-1-0. 77-l-O.OH-0. 08=8.86 sec 

k = F520 0.35 sec 1.14 sec 28.4-1-0.04=28.44 sec 39.3-1-9.1-1-0.05-1-0.49=48.94 sec 

Even if the theoretical complexity is linear in the size of k it is clear that, in 
practice, the behavior of the algorithm is not linear in log(fc). Moreover, when we 
increase the size of fc, step 5 becomes the most consuming part of our algorithm. 
Hence, even if the new algorithm is efficient enough to solve the problem for a 
small base field fc, the problems become intractable when #fc > 5^°°. In the next 
paragraph we propose several optimizations to overcome this limitation. 



Optimizations The idea is to apply recursively the algorithm of section 
to perform the step 5: we split again the first of variable into two parts: X 
X'UY' = X'U [X42,a;2i,a;5i,a:i2]. 



Algo 16. 3I Original Step 5 



Recursive Step 5 



8.0-^0.77+0.01-1-0.08=8.86 sec 0.05-h0.41-h0.33-h0.01=0.8 sec 
39.3+9.1+0.05-^0.49=48.94 sec 0.12+1.53+2.44+0.01+0.02=4.1 sec 

0.13+2.46+7.16+0.01+0.01=9.78 sec 



fc = F510 
fc = F520 
fc = F540 



When fc = F3fc we obtain in step 3 of the algorithm 16.31 the following lexico- 
graphical Grobner basis: 



' 9 (xoi) of degree 70 

hi {xo2, xoi) = 51 (xoi) {xl2 H ) and gi of degree 39 

h2 {xo2,xqa) = a;j]2 H 

• • • polynomials in variables X31, X32, a;o2, 2:01 



and thus we can split gi (xoi) into two factors: 

gi (xqi) = (a;oi + aif (xoi + 02)^ • • • (xqi + a^f 

^ ^^°^\ = xoi {xoi + (3if (xoi + ^2)^ • • • {xoi + ^3)^ 
91 (xoi) 

Hence the polynomial gi (xqi) can be efficiently factorized when fc is big. 



6.6 Experimental results 



Programming language Workstation 

The experimental results have been obtained with several Xeon bi-processor 3.2 
Ghz, with 16 Gb of Ram. The instances of our problem have been generated us- 
ing the Magma software. We used the Magma version 2.14 for our computations. 
The F5 |Fau02j algorithm have been implemented in C language in the FGb soft- 
ware and we used this implementation for computing the first Grobncr base. All 
the other computations are performed under Magma including factorizing some 
univariate polynomials and computing Grobncr bases using the F4 algorithm. 

Table Notation 

The following notations are used in the tables of Fig.l and Fig. 2 below: 



— fc is the ground field, fc' D fc is the field extension. The practical behavior of 
our algorithm is strongly depending on the size of fc'; hence, since k is fixed, 
the practical depends strongly on the degree of the field extension [k' : fc]. 
In order to obtain consistent data in the following tables we keep only the 
case [k' : k] = 2. 

— T is the total CPU time (in seconds) for the whole algorithm. 

— Tgcii is the time for generating the Riemann equations and computing a 
valid level 2 theta null point (Magma). 

~ ^Grob is the sum of the Grobner bases computations (FGb and Magma). 

— Tpact is the sum of the Factorization steps (Magma) . 

— Ti is the total time of the algorithm excluding generating the equations: 
Ti=T- Tgc„. 



k 


k' 


Tgcii 




Tpact 


Ti 


T 




5100 


1.9 


2.7 


9.3 


12 


14 


570 


5140 


3.4 


3.3 


16.0 


19 


23 


5100 


5200 


19.5 


15.9 


116.7 


133 


152 


5150 


5300 


27.9 


16.8 


159.7 


177 


205 


5200 


5400 


141.3 


57.3 


401.0 


459 


600 


g250 


5500 


178.4 


62.1 


651.8 


715 


893 


g300 


56OO 


227.8 


86.7 


935.3 


1023 


1251 


g350 


5700 


674.8 


108.5 


1306.1 


1416 


2091 


5400 


58OO 


764.1 


100.5 


2411.3 


2513 


3277 


5450 


5900 


1144.0 165.3 


2451.3 


2619 


3763 


g500 


glOOO 


1070.1 185.4 


2990.0 


3177 


4247 


56OO 


51200 


1979.5 273.5 


4888.6 


5164 


7144 


5700 


51400 


3278.0 422.5 


6872.2 


7297 


10575 



Fig 1: Algorithm ^ = 3, characteristic of k is 5. 



k 


fc' 




Tbrob 


Tpact 


Ti 


T 


38U 


3itio 


3.6 


2.0 


0.4 


3 


7 


380 


3I60 


3.6 


2.0 


0.2 


3 


6 


3200 


3400 


29.0 


11.1 


6.9 


20 


49 


3600 


31200 


239.2 


36.2 


44.5 


88 


327 


3800 


3I600 


403.7 


50.6 


89.6 


150 


554 


31000 


32000 


591.8 


61.8 


151.0 


225 


816 


31500 


33000 


2122.0 


137.7 


474.5 


666 


2788 


33000 


36000 


11219.9 396.3 3229.6 3704 


14923 



Fig 2: Algorithm I — '6, characteristic of k is 3. 



Interpretation of the results 

— In characteristic 3, the hardest part is the generation of the equations and the 
computation of a vahd level 2 theta null point: Tccn ~ T. In characteristic, 
5 we have T « 3 Tcen. 

— The most consuming part in algorithm described in 16.31 is the univariate 
factorization. Moreover due to the implementation in Magma Tpact is not 
really linear in the size of k. 

— The algorithm is much more efficient in characteristic 3 since: 

• All the solutions occur with some multiplicity, hence we have to deal 
with non-square-free polynomials. As a consequence, the degree of the 
univariate polynomials can be decreased by taking the square-free part 
of the polynomials. 

• The corresponding Grobner bases are in not in shape-position: as explain 
in section [6.41 we can split the univariate polynomial by taking a gcd. 

— The algorithm is very efficient since we can completely find the solutions 
of the ideal J for sizes of the base field k = 2>^^^^ 01 k = 5^*^° which are 
interesting for point counting application. 

7 Conclusion 

In this paper, we have described an algorithm to compute modular correspon- 
dences in the coordinate system provided by the theta null points of abelian 
varieties together with a theta structure. As an application, this algorithm can 
be used to speed up the initialisation phase of a point counting algorithm |CL08| . 
The main part of the algorithm is the resolution of an algebraic system for 
which we have designed a specific Grobner basis algorithm. Our algorithm takes 
advantage of the structure of the algebraic system in order to speed up the 
resolution. We remark that this special structure comes from the action of the 
automorphisms of the theta group on the solutions of the system which has a 
nice geometric interpretation. In particular we were able count the solutions of 
the system and to identify which one correspond to valid theta null points. 
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